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The young star cluster system of the Antennae galaxies 
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Abstract 

The study of young star cluster (YSC) systems, pref- 
erentially in starburst and merging galaxies, has seen 
great interest in the recent past, as it provides im- 
portant input to models of star formation. However, 
even some basic properties (like the luminosity function 
[LF]) of YSC systems are still under debate. Here we 
study the photometric properties of the YSC system in 
the nearest major merger system, the Antennae galax- 
ies. We find evidence for the existence of a statistically 
significant turnover in the LF. 
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1 Introduction 

Star clusters (SCs) form nearly instantaneously through 
the collapse of giant molecular gas clouds. Hence, 
all stars within a SC are approximately coeval, share 
the same chemical composition, and therefore rep- 
resent a simple stellar population. A small num- 
ber of parameters, in particular their initial chemi- 
cal composition and initial stellar mass function, are 
enough to describe their colour and luminosity evolu- 
tion on the basis of a given se t of stellar evolution- 
ary tracks or isochrones (e.g. iLeitherer et alJ [1999; 



Peter Anders 
Uta Fritze 
Richard de Grijs 

1 Sterrenkundig Instituut, Universiteit Utrecht, P.O. Box 80000, 
NL-3508 TA Utrecht, The Netherlands 

2 Centre for Astrophysics Research, University of Hertfordshire, 
College Lane, Hatfield ALIO 9AB, UK 

3 Departmcnt of Physics & Astronomy, The University of 
Sheffield, Hicks Building, Hounsfield Road, Sheffield S3 7RH, UK 



Anders and Fritze-v. Alvenslebenll2003 ; Bruzual and Char lotJ 



2003). Therefore, observed spectrophotometric proper- 
ties of SCs are relatively easy and straightforward to 
interpret. 

SC formation is a major mode of all star forma- 
tion, and possibly even the dominant mode in strong 
starbursts triggered in gas-rich gal axy mergers (e.g., 



Meurerl Il995t Ide Grijs et al.M200 3a). In addition, as 



SCs inherit and conserve the chemical abundances at 
the place and time of their birth up to old ages, they 
are excellent tracers of their parent galaxy's properties 
in terms of star formation and chemical enrichment his- 
tory. 

One of the most basic and commonly used diagnos- 
tics to explore the properties of entire SC systems is 
their LF. While for old globular cluster systems the 
Gaussian shape of their LFs (and indeed mass func- 
tions ) is well esta blished (see e.g. lAshman and Zep3 
ll99Rt iHarrisl 199ll ). the situation for young SCs is still 
under discussion. While mainly LFs consistent with 
a power-law, resembling the MF of nearby molecu- 
lar cloud cores, are quoted for YSC systems ranging 
from Galactic open cluster to YSCs in major mergers 



g. van den Bergh and Lafontainell984 ; Hunter et al 



l2003HSchweizer and Seitzerll998HWhitmore et al.lll999l | 

see lde Grijs et al.ll2003bl for a recent compilation) , some 

studies find deviations from a power-law or direct evi- 

i 1 

dence for Gaussia n distributions (Ide Griis and Anders 
l2006t iFritze-v. AlvensleberJl999t iGoudfrooij et alJ20 04) J 

To contribute to our understanding of this astro- 
physical question, which ties directly to the fundamen- 
tal physical conditions of star formation in general, we 
performed an analysis of the young (~ — 100 Myr) 
SC system formed during the, still ongoing, merging 
process in the nearest major merger of two giant gas- 
rich spiral galaxies, the Antennae galaxies, i.e. NGC 
4038/39. 
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2 Photometry 

We reanalysed the most homogeneous broad-band 
dataset including the U band of the Antennae galaxies 
available (providing photometry in UBVI), obtained 
using HST/WFPC2 as part of programme GO-5962 
(PI B. Whitmore). 

We select our clusters based on: 

• a flux threshold compared to its surrounding 

• cross-correlation between filters (i.e. a cluster must 
be detected in several bands) 

• photometric errors < 0.2 mag in ea ch band 

• size measurements using BAOlab (|Larsenlll999h 



1 cr 



Selecting only clusters in the (PSF-corrected) size 
range R\/2 — 5 to 25 pc (to avoid contamination by 
single stars and cluster complexes, respectively) yields 
a sample of 365 clusters, satisfying all selection criteria. 
Utilising the cluster sizes we app ly size-depend e nt ap er- 
ture corrections, as described in Anders et al. ( 2006I ). 
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Fig. 1 Completeness fractions for artificial clusters with 
FWHM = ~ 10 pc for a variety of selection criteria (see 
text for details) 



3 Completeness 

In this study we are primarily interested in the LF of 
the young star clusters. A consistent and reliable de- 
termination of the completeness fraction as function of 
cluster magnitude and all other relevant selection cri- 
teria/effects is therefore of prime importance. 

We investigated the completeness for a number of 
restrictions, taking successively more cluster selection 
criteria into account. Some of these results for artificial 
clusters with FWHM = ~ 10 pc in a cluster-rich region 
in NGC 4038 are shown in Fig. E 

First, we determined the completeness in each band 
independently (labelled "U","B", "V" and "I" in Fig. 
HJ. Taking the cross-correlation into account (labelled 
"XID") reduced the completeness to slightly below the 
most limiting single-band completeness. By apply- 
ing the photometric uncertainties restriction (labelled 
"XID + photo") the 50% completeness magnitude is 
decreased by about ~ 0.4 mag. 

Using stars instead of extended objects as test 
sources artificially increases the 50% completeness mag- 
nitude by ~ 1 mag. 

Computing completeness functions for 2 regions 
(characterised by different background and cluster den- 
sity levels) and 2 sizes (characteristic for our observed 
cluster sample) and imposing the size restriction onto 
the artificial clusters (labelled "XID photo size") al- 
lowed us to attribute realistic completeness fractions to 
each individual cluster, strengthening the subsequent 
analysis of the LFs. Imposing the size restrictions 
decrease the completeness by up to ~ 1 mag, while 



different cluster sizes and regions within the galaxies 
contribute completeness changes of ~ 0.6 and ~ 1 mag, 
respectively. This very complex dependence of the com- 
pleteness fractions on a number of input parameters 
clearly shows the necessity to model the complete- 
ness functions as realistically as possible, taking 
into account all cluster selection criteria. For the 
same region as in Fig. [TJ results for different cluster 
sizes and selection criteria are shown in Fig. [2] 

For add itional details of the completeness determi- 
nation, see 



Anders et all (120071 ). 



4 The model to analyse the luminosity 
functions 

In collaboration with statisticians from the University 
of Gottingen we developed a suite of statistical tools. 
This allows to fit the LFs of our young star clusters with 
either a Gaussian distribution or a power-law, access 
the accuracy of the fits and the fitted parameters, and 
establish the statistical significance of differences in the 
goodness-of-fit between the two fitted test distributions. 

The fitting of a observed LF takes the photomet- 
ric error and appropriate completeness fraction of each 
individual cluster into account, and independently fits 
either a Gaussian distribution or a power-law to the 
LF, yielding a set of parameters characterising the best- 
fitting model and a likelihood parameter for either test 
distribution. 

We find a strong superiority (measured by the ratio 
of the likelihoods for Gaussian and power-law fits) of the 



YSCs in the Antennae galaxies 
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1 pix XID phot 
- 1 pix XID phot size 

• 2pix XID 
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Fig. 2 Completeness fractions for artificial clusters with 
FWHM = ~ 10 pc (= lpixel) and FWHM = ~ 20 pc (= 
2pixel) for a variety of selection criteria (see text for details) 



Gaussian fit compared to the power-law fit (see Sect. E|) . 
To validate this result, we have to test its statistical sig- 
nificance, i.e. to determine the probability to achieve 
such strong superiority of the Gaussian distribution if 
the true underlying LF is a power-law, distorted by the 
error distribution and the completeness function. In a 
Monte-Carlo approach we draw artificial observations 
from the best-fitting power-law, use our statistics tools 
to fit them again with Gaussian and power-law distribu- 
tion, and determine the likelihood ratio for these artifi- 
cial observations. The likelihood ratio of the observed 
distribution compared to the distribution of likelihood 
ratios from our artificial tests then yields the probabil- 
ity that the determined superiority of the Gaussian fits 
is consistent with an underlying power-law distribution. 

The uncertainties of the best-fitting model parame- 
ters are estimated by bootstrapping. 

For additi o nal d etails of the statistical model, see 
Anders et all (|2007h . 



5 Results 

As example, the LF in the V band for our YSC sam- 
ple in the Antennae galaxies is shown in Fig. From 
visual inspection, the best fit with an underlying Gaus- 
sian distribution appears to represent the data better 
than the best fit with an underlying power-law distri- 
bution. 

To quantify this, we perform our likelihood ratio test 
and the Monte-Carlo analysis. For the observed distri- 
butions we find a likelihood ratio value of 23.1 (where 
larger values are equivalent to a stronger superiority of 




Fig. 3 V-band luminosity functions for our Antennae clus- 
ters and fit: Observed (green solid line), best fit (black solid) 
and underlying distribution of best fit (black dashed). Fit- 
ted test distributions are: Gaussian (left panel) and power- 
law (right panel) 



the Gaussian fit when compared to the power-law fit). 
The Monte-Carlo analysis with 1000 test realisations 
drawn from the best-fitting power-law distribution re- 
sults in a maximum likelihood ratio of 11.4, hence none 
of the test distributions can reproduce a superiority of 
the Gaussian distribution as strong as observed. Util- 
ising the properties of Bernoulli-distributed variables, 
this corresponds to a probability < 0.5% that the un- 
derlying distribution is still consistent with a power- 
law. 

This result is valid also for the other bands (except 
the I band, for which the observations are significantly 
shallower), and for several age and siz e subsamples. For 
more details, see lAnders et al.l ([20071 ). 
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